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Abstract 

The Neurachninae is the only grass lineage known to contain C 3 , C 4 , and C 3 -C 4 intermediate species, and as such 
has been suggested as a model system for studies of photosynthetic pathway evolution in the Poaceae; however, a 
lack of a robust phylogenetic framework has hindered this possibility. In this study, plastid and nuclear markers were 
used to reconstruct evolutionary relationships among Neurachninae species. In addition, photosynthetic types were 
determined with carbon isotope ratios, and genome sizes with flow cytometry. A high frequency of autopolyploidy 
was found in the Neurachninae, including in Neurachne munroi F.Muell. and Paraneurachne muelleri S.T.Blake, which 
independently evolved C 4 photosynthesis. Phylogenetic analyses also showed that following their separate C 4 origins, 
these two taxa exchanged a gene encoding the C 4 form of phosphoeno/pyruvate carboxylase. The C 3 -C 4 intermedi- 
ate Neurachne minor S.T.Blake is phylogenetically distinct from the two C 4 lineages, indicating that intermediacy in 
this species evolved separately from transitional stages preceding C 4 origins. The Neurachninae shows a substantial 
capacity to evolve new photosynthetic pathways repeatedly. Enablers of these transitions might include anatomical 
pre-conditions in the C 3 ancestor, and frequent autopolyploidization. Transfer of key C 4 genetic elements between 
independently evolved C 4 taxa may have also facilitated a rapid adaptation of photosynthesis in these grasses that 
had to survive in the harsh climate appearing during the late Pliocene in Australia. 

Keywords: C 4 grass evolution, C 4 photosynthesis, C 3 -C 4 intermediate, grass phytogeny, lateral gene transfer, Neurachne, 
Paraneurachne, polyploidy 



Introduction 

Despite its complexity, the C 4 photosynthetic pathway has C 4 lineages have been postulated in the PACMAD clade (Grass 

evolved independently >62 times in flowering plants (Sage et ah, Phylogeny Working Group II, 2012). Such a clustering of C 4 ori- 

2011), thus constituting a striking example of convergent evolu- gins is also observed in other groups, with six independent lin- 

tion. It is especially prevalent in grasses, where 22-24 distinct eages in the sedges (Cyperaceae) and 23 in the Caryophyllales 
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(Sage et al., 201 1). These patterns indicate that some plant groups 
have a higher propensity for C 4 photosynthesis evolution, which 
may reflect ecological, genomic, and/or anatomical factors that 
facilitate the acquisition of novel traits (Sage, 200 1 ; Marshall 
et al, 2007; McKown and Dengler, 2007; Christin et al, 2011; 
Edwards and Ogburn, 2012). Leading environmental factors pro- 
moting C 4 evolution are low atmospheric C0 2 , heat, drought and 
salinity, often in combination (Sage et al., 2012). Anatomical 
factors include high vein density, which may be common in dry 
environments and certain taxonomic groups such as the grasses 
(Ehleringer et al, 1997; Sage et al, 2012). 

The evolution of C 4 photosynthesis is best studied using closely 
related taxa with different photosynthetic types, as in the eud- 
icot groups Flaveria, Cleome, Molluginaceae, and Heliotropium 
(McKown et ah, 2005; Marshall et al., 2007; Feodorova et al., 
2010; Christin et al., 2011; Muhaidat et al., 2011). These groups, 
however, have limited utility for understanding the origins of the 
pathway in the grasses, where half of all C 4 species occur. C 4 
grasses are the most successful group of C 4 plants on Earth, dom- 
inating the 23% of global primary productivity attributable to C 4 
vegetation, and comprising the vast majority of C 4 plants in agri- 
cultural use (Brown, 1999; Still et al., 2003). Efforts to engineer 
the C 4 pathway into C 3 crops to take advantage of the superior 
productivity of C 4 photosynthesis have been directed towards 
the grasses rice and wheat (http://irri.org/c4rice); however, these 
efforts are hindered by the lack of a model group for studying C 4 
evolution in the Poaceae. Such a model group could be exploited 
to identify the genetic changes that occurred during C 4 evolu- 
tion, as well as elucidating the order in which the individual 
traits of the pathway were assembled. The major C 4 lineages of 
grasses (e.g. Chloridoideae and Andropogoneae) are composed 
of numerous and ecologically successful C 4 species, but lack 
C 3 -C 4 intermediate species, and are only distantly related to C 3 
taxa (Christin et al., 2008; Grass Phylogeny Working Group II, 
2012). Photosynthetic variation exists in some small groups of 
grasses, notably Steinchisma, which contains C 3 -C 4 intermedi- 
ates and C 3 species, and Alloteropsis , which has both C 3 and C 4 
taxa (Duvall et al., 2003; Ibrahim et al., 2009; Christin et al., 
2010), but only one grass clade has been identified that contains 
C 3 , C 4 , and C 3 -C 4 intermediate species. This is the Neurachne/ 
Thyridolepis clade (in the subtribe Neurachninae; Morrone 
et al., 2012), a group of three genera and 11 species endemic to 
Australia (Blake, 1972; Macfarlane, 2007). 

Within the Neurachninae, the genus Neurachne includes one 
C 4 species (/V. munroi), five C 3 species, and one C 3 -C 4 interme- 
diate species, N. minor (Hattersley et al., 1982, 1986; Hattersley 
and Roksandic, 1983; Macfarlane, 2007). The monospecific 
genus Paranenrachne (P. muelleri) is C 4 , while the other genera 
(Thyridolepsis, Ancistrachne, Cleistochloa, and Calyptochloa) 
are C 3 . The clade belongs to the Panicoideae subfamily, which 
encompasses the vast majority of C 4 grass lineages, and thus 
appears especially prone to transitions from C 3 to C 4 biochemistry 
(Grass Phylogeny Working Group II, 2012). Extensive work in the 
1980s characterized the anatomy, biochemistry, and physiology of 
Neurachninae species (Hattersley et al., 1982, 1986; Hattersley and 
Roksandic, 1983; Hattersley and Stone, 1986; Moore andEdwards, 
1989), indicating high potential for this group to serve as a model 
for C 4 grass evolution. However, to make evolutionary inferences, 



it is necessary to have a well-resolved, species-level phylogeny. 
Such a phylogeny was not available for the Neurachninae, as 
only some members of the group have been analysed with a small 
number of molecular markers (Hudson et al., 1990; Christin et al., 
2008; Grass Phylogeny Working Group II, 2012). 

The objective of the present study was a reconstruction of 
the evolutionary history of the Neurachninae, with an emphasis 
on photosynthetic pathway evolution. Multiple accessions per 
species were sampled, and phylogenetic analyses of plastid as 
well as nuclear markers, photosynthetic pathway identification, 
and measures of genome size were included. This comparative 
approach established the phylogenetic relationships between 
Neurachninae species, and revealed the genome dynamics of the 
group. The outcomes also highlight the diversity of photosyn- 
thetic transitions in the Neurachninae, and provide a solid foun- 
dation for future studies aimed at elucidating the anatomical and 
molecular mechanisms underlying these transitions. 

Materials and methods 

Plant material 

Live, field-collected individuals of Neurachninae species were sampled 
for genome size analyses using flow cytometry (Supplementary Table 
S 1 available at JXB online), while carbon isotope ratios were determined 
using multiple herbarium samples of each species (Supplementary 
Appendix SI). Multiple individuals of each Neurachninae species were 
also sampled for phylogenetic analyses; these were herbarium speci- 
mens or plants collected from the field (Supplementary Tables SI, S2). 

Carbon isotope ratios (8 13 C) 

Carbon isotope values were determined for all described Neurachne, 
Paraneurachne, and Thyridolepis species (Supplementary Appendix SI 
at JXB online) by the University of California, Davis Stable Isotope 
Facility (http://stableisotopefacility.ucdavis.edu). Some of the samples 
were previously assayed for 6 13 C (Hattersley and Roksandic, 1983; 
Hattersley et al., 1986). 

Determination of genome size 

Nuclei were simultaneously released from fresh leaf material from a 
Neurachninae species and a calibration standard [Raphanus sativus cv. 
Saxa (2C DNAcontent= 1.11 pg; Dolezel et al., 1992) or Lepidosperma 
gibsonii (2C DNA content = 0.56 pg; M. Wallace, unpublished results)] 
by chopping with a razor blade in cold buffer (Roberts, 2007). Samples 
were filtered, and the nuclei were stained with propidium iodide 
(Roberts, 2007) and analysed using a BD FACSCanto II flow cytom- 
eter (Becton Dickinson, Franklin Lakes, NJ, USA) with 488 nm excita- 
tion and 585/424 band pass filter detection. Data were collected with 
BD FACSDIVA software (v. 5.0.2) and analysed using FlowJo v.7.6.3 
software (Tree Star Inc., Ashland, OR, USA). At least 5000 nuclei per 
day were analysed from each individual on three separate days. The 
genome size of an individual was then calculated as the average of these 
three estimates. Genome size and standard deviation (SD) of conspecific 
homoploids (Supplementary Table SI at JXB online) were calculated 
using the average genome sizes of the individual homoploids. Estimates 
were only included if the peak heights of the sample and standard were 
similar, and the coefficient of variation was <5%. Conversion of pg 
DNA to Mbp DNA followed Dolezel et al. (2003). 

Amplification of plastid markers and ITS 

Five plastid markers (ndhF, trnK/matK, rpoC2, rpll6, and trnLF) that 
have been used previously to investigate relationships among grasses 
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(Duvall et al., 2003; Salariato et al., 2010; Grass Phylogeny Working 
Group H, 2012) were amplified from genomic DNA (gDNA) extracted 
from herbarium samples (MP FastDNA SPIN Kit and FastPrep 
Instrument; MP Biomedicals LLC, Solon, OH, USA). For ndhF and 
trnK/matK, previously published primers were used (Grass Phylogeny 
Working Group II, 2012), while primers for the three other mark- 
ers were designed in the conserved regions of Panicoideae sequences 
downloaded from GenBank (Supplementary Table S3 at JXB online). 
The nuclear region encompassing the internal transcribed spacer 1, 5.8S 
rRNA, and internal transcribed spacer 2 (ITS) was amplified with the 
universal primers ITS4_rev and ITS5_for (Supplementary Table S3; 
White et al., 1990). However, in half the samples, these primers ampli- 
fied endophytic fungal genes. A new forward primer, specific to grasses 
(ITSgrassesfor; Supplementary Table S3), was then used in combina- 
tion with the ITS4 reverse primer. The plastid and ITS markers were 
amplified in overlapping fragments of -300-700 bp, with the protocol 
described in Grass Phylogeny Working Group II (2012). Single nucleo- 
tide polymorphisms were detected in most of the ITS sequences. These 
were coded as ambiguous characters in the phylogenetic analyses. All 
sequences have been submitted to GenBank (Supplementary Table S2). 

Amplification of low-copy nuclear markers 

As the gDNA extracted from herbarium specimens was too degraded 
to amplify low-copy nuclear markers, gDNA was isolated (DNeasy 
Plant Mini Kit, QIAGEN, Hilden, Germany) from selected accessions 
for which fresh material was available (Supplementary Table S4 at 
JXB online). Fragments of the genes encoding waxy and arodeh were 
amplified with previously designed primers (Christin et al., 20126). 
The ppc-B2 lineage of phosphoOTo/pyruvate carboxylase (PEPC), 
which contains most C 4 -specific forms in grasses as well as non-C 4 
orthologues (Christin et al., 2007), was isolated in one fragment or two 
overlapping segments, using previously published primers (Christin 
et al., 2007, 20126). These three markers were isolated as described 
previously (Christin et al., 20126), but with an annealing temperature of 
5 1 °C and an extension time of 2 min. All sequences have been submit- 
ted to GenBank (Supplementary Table S4). 

Phylogenetic analyses 

The veracity of the sequence data was monitored throughout the analy- 
ses. The congruence between accessions of the same species and the 
lack of identical sequences for any of the plastid or nuclear markers, 
including the ppc-B2 genes, from different Neurachninae species indi- 
cated no cross-contamination or misidentification of the samples. 

Sequences were aligned with MUSCLE (Edgar, 2004), and the align- 
ments were manually refined. Outgroups were included based on pre- 
vious work (Christin et al., 20126; Grass Phylogeny Working Group 
II, 2012). Best substitution models were determined by hierarchical 
likelihood ratio tests. In all cases, the general time reversible substi- 
tution model with a gamma shape parameter (GTR+G) was selected. 
Phylogenetic trees were obtained by Bayesian inference with MrBayes 
3.2 (Ronquist and Huelsenbeck, 2003). Two analyses, each consisting of 
four parallel chains, were run for 7 000 000 generations after a burn-in 
period of 3 000 000, sampling a tree each 1000 generations. A consen- 
sus tree was computed on the 14 000 sampled trees. Phylogenetic trees 
were also inferred under a maximum likelihood (ML) criterion with the 
software PhyML (Guindon and Gascuel, 2003). Branch support of the 
ML trees was evaluated with 100 bootstrap pseudo-replications. 

Phylogenetic trees were inferred for concatenated plastid data and ITS 
separately, and a further estimate of phylogeny was carried out by con- 
catenating these two data sets. The three low-copy nuclear markers were 
analysed separately. For waxy, parts of the sequences that consisted of 
simple sequence repeats or short tandem repeats were difficult to align 
and were removed from the data set, resulting in the elimination of 52 
aligned base pairs. To decrease the possibility of a bias due to conver- 
gent evolution (Christin et al., 2007, 20 12a), ppc-B2 phylogeny was also 
inferred from introns only. The significance of topological incongruence 
between the different markers was evaluated by Shimodaira-Hasegawa 



(S-H) tests with multiple-comparison correction (Shimodaira and 
Hasegawa, 1999), as implemented in baseml software (Yang, 2007). 

Molecular dating 

The divergence times of species within the Neurachninae were esti- 
mated with a molecular dating approach, which allows for rate variation 
among branches, following the recommendations from Rutschmann 
(2006), as described in Christin et al. (2008). Molecular dating requires 
time constraints set a priori on some nodes of the phylogeny as mini- 
mum or maximum ages. In the absence of a reliable fossil record for 
this group, the root of the Neurachninae was fixed to 1 1 Million years 
ago (Ma), a preliminary estimate obtained with a larger phylogeny, 
and without taking phytolith fossils into account (P. -A. Christin et al., 
unpublished). If phytoliths are taken into account, this node would prob- 
ably be moved to ~20 Ma (P.-A.Christin et al., unpublished); however, 
as the same constraint was set for all markers, the relative ages can be 
compared among the different markers. 

The topology inferred from combined plastid markers and ITS was 
used for molecular dating as it was better resolved and more complete 
than those of individual markers. However, the estimation of ages was 
based only on plastid markers to avoid strong variation in evolutionary 
rates due to differences between plastid and nuclear markers. Additional 
molecular dating analyses were carried out from the ITS, as well as 
each of the low-copy nuclear genes, with topologies compatible with 
the inferred species relationships. For these markers, an outgroup from 
outside the Neurachninae was added to the data set and used to estimate 
the first divergence time of the Neurachninae, but was removed during 
the molecular dating analyses. 

Results 

Photosynthetic pathway variation in the Neurachninae 

Leaf 6 I3 C assays confirmed previous photosynthetic pathway 
determinations in the Neurachninae (Hattersley and Roksandic, 
1983; Hattersley et al., 1986; Macfarlane, 2007). Mean 6 I3 C val- 
ues near — 13%o in A', munroi and P. muelleri demonstrate that 
these species are C 4 (Table 1). All other Neurachninae members, 
including the C 3 -C 4 intermediate A', minor, showed 8 13 C val- 
ues between -24%o and -29%o, which are typical of C 3 plants 
(Table 1). 

Neurachninae genome size analyses 

The 2C DNA estimates varied 4-fold from 0.98 ±0.01 pg DNA 
for the diploid N. lanigem S.T.Blake to 3.97±0.04 pg DNA for 
the tetraploid P. muelleri (Table 1). These values are well within 
the range observed in the Poaceae (Bennett and Leitch, 2010, 
Plant DNAC-values database, http://data.kew.org/cvalues/), and 
agree with previous estimates from the Neurachninae based on 
2C values or chromosome counts (Prendergast and Hattersley, 
1985; Morgan and Westoby, 2005; Macfarlane, 2007). 

Intraspecific ploidy variation was found in 7Y. alopecuroidea 
R.Br. (4x, 5x, 6x, and 7x), and although multiple ploidy levels 
have been observed previously for N. lanigem and N. munroi 
(Prendergast and Hattersley, 1985), only diploid and tetraploid 
forms of these species, respectively, were detected in this study. 

Phylogenetic relationships among the Neurachninae 

All plastid and nuclear markers showed that Neurachne and 
Paraneurachne species formed a strongly supported monophyletic 
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Table 1. Carbon isotope ratios (b 13 C) and genome size data for Neurachninae species. Carbon isotope ratios are average values 
±standard deviation (SD) determined from multiple herbarium samples of the same species (Supplementary Appendix S1 at JXB online). 
Photosynthetic pathway designations are shown based on these results and Hattersley ef al. (1 986) for the C 3 -C 4 species N. minor. 
Genome size measurements ±SD were calculated from the average genome sizes of individual homoploids, which were measured from 
fresh leaf material on three consecutive days (see the Materials and Methods for further detail). 



Species 


6 13 C ±SD 


2C ±SD (pg) 


2C (Mbp) 


1Cx ±SD (pg) 


2n 


Ploidy 


Neurachne alopecuroidea (C 3 ) 


-27.1 +1.5 


2.24 + 0.02 


2190 


0.56 + 0.01 


36 


Ax 






O ~7~7 f~\ f~\r? 

2.77 + 0.05 


271 0 


0.55 + 0.01 


45 a 


5x H 






3.33 + 0.04 


3260 


0.55 + 0.01 


54 


6x 






3.81 


3730 


0.54 


63 a 


Ix 1 ' 


N. annularis (C 3 ) 


-24.8 + 0 


1 .30 + 0.01 


1270 


0.65 + 0.01 


18 


2x 


N. lanigera (C 3 ) 


-28.4 + 1.2 


0.98 + 0.01 


958 


0.49 + 0.01 


18 


2x 






ND 


ND 


ND 


36 


Ax 


N. minor (C 3 -C 4 ) 


-27.1+1.0 


2.76 


2700 


0.69 


36 


Ax 


N. munroi (C 4 ) 


-12.7 + 1.0 


ND 


ND 


ND 


18 


2x 






3.61 


3530 


0.60 


36 


Ax 






ND 


ND 


ND 


54 


6x 


N. queenslandica (C 3 ) 


-27.1 +0.4 


ND 


ND 


ND 


54 


6x 


N. tenuifolia (C 3 ) 


-25.7 + 0.5 


1.48 


1450 


0.74 


18 


2x 


Paraneurachne muelleri (C 4 ) 


-13.3 + 0.4 


3.97 + 0.03 


3880 


0.99 + 0.01 


36 


Ax 


Thyridolepis mitchelliana (C 3 ) 


-26.6 + 1.4 


1.32 + 0.01 


1290 


0.66 + 0.01 


18 


2x 


T. multiculmis (C 3 ) 


-26.3 + 1.8 


2.50 + 0.02 


2450 


0.63 + 0.01 


36 


Ax 


T. xerophila (C 3 ) 


-27.6 + 0.6 


ND 


ND 


ND 


18 


2x 



ND, not determined. 

a Ploidy and chromosome numbers are inferred from the flow cytometry results of this study; other values are from Prendergast and Hattersley 
(1985) and Macfarlane (2007). 



queenslandica, therefore, may be a natural hybrid between C 3 
and C 4 parents, although further support for this hypothesis is 
required as C 4 photosynthesis may have evolved in P. muelleri 
after the hybridization event. 

The phylogeny inferred from combined plastid and ITS mark- 
ers (plastid+ITS) was strongly resolved (Fig. 1), and was thus 
considered representative of the species relationships. In this 
phylogeny, the monophyly of all species was strongly supported, 
with the exception of Thyridolepis taxa. 

The relationships between species inferred from arodeh 
(Supplementary Fig. S3 at JXB online) and waxy (Supplementary 
Fig. S4) were poorly supported, but nevertheless compatible 
with the plastid+ITS phylogeny. Forcing the topologies to that 
deduced from plastid+ITS did not significantly decrease the 
likelihood (S-H test, difference of log-likelihoods = -2.786 
and -0.822, respectively, P=0.257 and 0.269, respectively). 
All arodeh and waxy sequences isolated from the same species 
were monophyletic, with the exception of the sequences encod- 
ing these markers isolated from T. multiculmis S.T.Blake. Two 
sequences encoding both waxy and arodeh were isolated from 
this species, and in both cases one sequence was more similar 
to that from T. mitchelliana S.T.Blake than it was to the other 
sequence isolated from T. multiculmis. 

Evolutionary history of ppc-B2 in the Neurachninae 

The coding sequences encompassing exons 5-10 (1786bp; 270 
parsimoniously informative sites within the Neurachninae) of 
the Neurachninae ppc-B2 genes were placed in a larger data set 



clade, which was sister to Thyridolepis spp., and the monospecific 
Paraneurachne was nested inside Neurachne (Supplementary 
Figs S1-S4 at JXB online). The original name of this species, 
Neurachne muelleri Hack. (Hackel, 1895), should thus be 
resurrected. 

Based on plastid markers, each species of the Neurachne/ 
Paraneurachne group was strongly supported as monophyletic, 
with the exception of A', lanigera, which was poorly resolved 
with respect to A', alopecuroidea (Supplementary Fig. S 1 at JXB 
online). The two C 4 taxa (N. munroi and P. muelleri) were not 
closely related, and the C 3 -C 4 intermediate A^. minors/as not sis- 
ter to either of the C 4 taxa (Supplementary Fig. SI). 

The topology of the phylogeny inferred from ITS 
(Supplementary Fig. S2 at JXB online) was almost identical to 
that from plastid markers. The only exception was A", queens- 
landica S.T.Blake, which was positioned as sister to A', alope- 
curoidea and not P. muelleri. Forcing A^. queenslandica to be 
sister to A', alopecuroidea in the phylogeny inferred from plas- 
tid markers led to a significant decrease of likelihood (S-H test, 
difference of log-likelihoods= -55.801, P < 0.001), as did forc- 
ing it to be sister to P. muelleri in the ITS phylogeny (S-H test, 
difference of log-likelihoods= -35.332, P < 0.005). This incon- 
gruence between the plastid and ITS inferred phylogenies thus 
appears real and is symptomatic of incomplete lineage sorting 
or reticulate evolution. Since N. queenslandica is hexaploid 
(Prendergast and Hattersley, 1985), an allopolyploid origin is 
likely, with the female parent related to the P. muelleri lineage 
and the male parent to the N. alopecuroidealN. lanigera line- 
age (assuming maternal chloroplast inheritance). Neurachne 
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encompassing many Panicoideae ppc-B2 sequences (Christin 
et al., 20126). In this inferred phylogeny, the Neurachninae 
sequences were positioned near those of affiliated C 3 Paniceae, 
and were monophyletic, with the exception of one Ancistrachne 
sequence, which was positioned outside the Neurachninae and 
was therefore excluded from further analyses. 

Inside the Neurachninae, despite a lack of resolution of the 
deeper nodes, the phylogenetic relationships deduced from 
ppc-B2 sequences were almost identical to those deduced from 
other markers (Supplementary Fig. S5 at JXB online). The only 
exceptions were the sequences from P. muelleri and N. mun- 
roi. Three groups of highly divergent ppc-B2 sequences were 
retrieved from P. muelleri. One of these groups was supported 
as sister to N. alopecuroidealN. lanigera ppc-B2 sequences, as 
expected based on other markers (Fig. 1). The exact positions 
of the other two clusters of P. muelleri ppc-B2 sequences, how- 
ever, were not resolved with confidence, although one of them 
is strongly supported as sister to sequences from N. munroi 
(Supplementary Fig. S5, subclade highlighted in red). The close 
relationship between N. munroi ppc-B2 sequences and some 
ppc-B2 sequences from P. muelleri was strongly supported even 
when a phylogeny was inferred from the five introns only (774 
aligned base pairs, including 129 parsimoniously informative 
sites within the Neurachninae; Supplementary Fig. S6, subclade 
highlighted in red). 

Different hypotheses reconciling ppc-B2 phylogeny with those 
inferred from other markers were evaluated using the sequence 



data (Fig. 2). Forcing the monophyly of all P. muelleri ppc-B2 
sequences and placing them as expected for the species, and 
placing the N. munroi sequences in the predicted species posi- 
tion — the hypothesis of gene transmission that followed species 
genealogy — very strongly decreased the likelihood of the data 
(S-H test, difference of log-likelihoods= -246.516, P < 0.0001). 
Thus the close relationship between some of the P. muelleri ppc- 
B2 sequences and the N. munroi ppc-B2 sequences is unambigu- 
ous. A scenario in which a gene transfer occurred from N. munroi 
to P. muelleri involved placing the closely related N. munroi and 
P. muelleri sequences as expected for N. munroi, and the other 
P. muelleri sequences as expected for this species (Fig. 2A). 
This treatment of the data did not significantly affect the likeli- 
hood (S-H test, difference of log-likelihoods= -5.542, P=0.621). 
The alternative hypothesis, a gene transfer from P. muelleri to 
N. munroi, was also examined, and involved moving all iV. mun- 
roi and P. muelleri ppc-B 2 sequences to the position expected for 
P. muelleri (Fig. 2B). Again, no significant effect on the likeli- 
hood was seen (S-H test, difference of log-likelihoods= -6.799, 
P=0.588). The hypothesis that the data set might encompass dif- 
ferent paralogues that appeared through gene duplication before 
the diversification of Neurachnel Paraneurachne was also tested. 
This implied placing P. muelleri sequences as sister to all other 
ParaneurachnelNeurachne sequences (Fig. 2C), which did not 
significantly decrease the likelihood (S-H test, difference of log- 
likelihoods= -9.391, P=0.485). This hypothesis, however, is not 
supported by the molecular dating analyses, which indicate that 
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Fig. 1. Phylogeny of the Neurachninae, photosynthetic types, and ploidy levels. The tree was obtained through Bayesian inference 
using the plastid markers ndhF, trnK/matK, trnLF, rpll6, and rpoC2, and the genomic region encoding the internal transcribed spacer 
1 , 5.8S rRNA, and internal transcribed spacer 2 (ITS). Bayesian support values and bootstrap values are indicated above and below 
branches, respectively. The tree was calibrated and branch lengths are proportional to divergence times, in million years ago (Ma). 
Photosynthetic types are indicated on the right, as are ploidy levels. Multiple branches to the same species denote different genes or 
alleles. 
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Fig. 2. Different evolutionary scenarios for ppc-B2, and divergence times within the Neurachninae. The inferred chronograms are 
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the divergence of the closely related N. munroi and P. muelleri 
ppc-B2 genes occurred long after the initial diversification of 
NeurachnelParaneurachne (Fig. 2C). Finally, it was considered 
whether the phylogenetic analyses might be biased by conver- 
gent amino acid changes in these two C 4 species (Christin et al., 
2007, 2012a). However, topological tests based only on introns 
yielded identical results, ruling out such a scenario. The data, 
therefore, show that at least some Neurachninae ppc-B2 genes 
were not transmitted following the species genealogy, but the 
direction of transmission, N. munroi to P. muelleri or vice versa, 
cannot be determined with confidence from the current data. 

Residues at positions putatively selected for C 4 -specific func- 
tion (Christin et al., 2007) were compared among the predicted 
amino acid sequences of the Neurachninae ppc-B2 sequences 



(Fig. 3). Of the 12 positions with the highest probability of hav- 
ing been positively selected in grasses (Christin et al., 2007), four 
were conserved in all members of the Neurachninae, with amino 
acid residues characteristic of non-C 4 ppc-B2 proteins (positions 
577, 637, 761, and 807). At the other positions, the predicted 
amino acid sequences encoded by the closely related N. munroi 
and P. muelleri ppc-B2 genes presented residues that character- 
ize independently evolved C 4 ppc-B2 genes (Fig. 3), with a ser- 
ine at position 780 that was shown to be a major determinant of 
C 4 -specific PEPC activity (Biasing et al., 2000; Svensson et al., 
2003). The other P. muelleri ppc-B2 sequences were predicted 
to encode some C 4 -characteristic residues (Fig. 3), but also con- 
tained several stop codons and deletions altering the reading 
frame, and are very probably non-functional pseudogenes. 
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Fig. 3. Evolutionary history of ppc-B2 in the Neurachninae. The calibrated phylogeny represents one of the hypotheses compatible with 
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Molecular dating of the Neurachninae 

The confidence intervals of estimated species divergence times 
based on plastid and nuclear markers overlapped (Fig. 4). Within 
most species, the divergence time between nuclear sequences 
overlaps with the crown age of the species. However, N. munroi 
arodeh sequences are estimated to have diverged long before the 
different N. munroi accessions. The occurrence of this pattern 
on a single gene can be explained by gene duplication or poly- 
ploidy followed by the loss of some gene copies. In addition, 
with multiple nuclear markers, the divergence of the different 
T. multiculmis sequences is estimated before the divergence of 
the different Thyridolepis species. This pattern, restricted to the 
tetraploid T. multiculmis, suggests an allopolyploid origin, with 
parents inside Thyridolepis, one of which was either not sam- 
pled or is extinct. The different hypotheses for the evolutionary 
history of ppc-B2 described above yielded similar age estimates 
(Fig- 2). 

Discussion 

Evolutionary origins of C 4 biochemistry in the 
Neurachninae 

The C 4 and C 3 -C 4 members of the Neurachninae are not sister 
species, but are separated in both plastid- and nuclear marker- 
based phylogenies by C 3 taxa (Fig. 1; Supplementary Figs 
S1-S4 at JXB online). In other groups containing species dem- 
onstrating different photosynthetic types, and for which robust 



phylogenies exist, some C 3 -C 4 species are directly sister to C 4 
taxa as in Flaveria (McKown et al., 2005), Cleome (Feodorova 
et al, 2010), Heliotropium (Hilger and Diane, 2003), and the 
Molluginaceae (Christin et ah, 201 1). This pattern suggests that 
the photosynthetic variation observed in extant taxa reflects a 
gradual transition from C 3 to C 3 -C 4 and finally to C 4 photosyn- 
thesis. However, several C 3 -C 4 taxa are not closely related to any 
C 4 species (Sage et al., 20 1 1 ), and many groups encompass a mix- 
ture of C 3 , C 3 -C 4 , and C 4 taxa, including the Neurachninae, but 
also other Flaveria (McKown et ah, 2005) and Molluginaceae 
(Christin et al., 2011) lineages, sedges (Roalson et al., 2010), and 
the Chenopodiaceae (Sage et al., 2011). 

The clustering of taxa using different photosynthetic pathways 
that is observed in the phylogenetic trees of the Neurachninae 
and other groups suggests a complex evolutionary history, with 
either multiple transitions from C 3 to C 3 -C 4 or C 4 photosyn- 
thesis, multiple losses of an ancestral C 4 or C 3 -C 4 type, or a 
combination of both (Christin et al., 2010). A simple scenario 
of a single C 4 origin followed by reversals to the C 3 (or C 3 -C 4 ) 
state can be ruled out for the Neurachninae by the analyses of 
genes encoding PEPC. During the evolution of this key C 4 path- 
way enzyme, mutations occurred in the genes that resulted in 
modifications to the catalytic and regulatory properties of the 
encoded proteins, allowing them to function optimally in C 4 
photosynthesis (Dong et al, 1998; Biasing et al., 2000). These 
mutations often occur at identical positions in the PEPC cod- 
ing sequences of distantly related species (Christin et al., 2007; 
Besnard et al, 2009). Several amino acid changes characteristic 
of C 4 -specific PEPC were predicted from the ppc-B2 genes of C 4 
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Neurachninae, including the key C 4 determinant alanine to serine 
mutation at position 780 (Fig. 3). If the C 4 -specific properties of 
PEPC appeared before the diversification of the Neurachninae, 
the ppc-B2 genes of non-C 4 taxa most probably would have 
maintained codons for some C 4 -characteristic amino acids. 
However, the ppc-B2 genes of the C 3 and C 3 -C 4 Neurachninae 
do not encode C 4 -specific amino acids, and are indistinguishable 
at these positions from species known never to have been in a C 4 
state (Fig. 3), strongly suggesting that the ppc-B2 genes of these 



non-C 4 taxa did not descend from genes encoding enzymes that 
were C 4 optimized. 

The C 4 -specific properties of the PEPC used by N. mun- 
roi and P. muelleri evolved after the divergence from non-C 4 
Neurachninae, but the C 4 PEPC genes of these two species 
strongly grouped together in the phylogeny, even when only 
introns were considered, arguing against genetic convergence as 
a phylogenetic bias (Supplementary Figs S5, S6 at JXB online). 
The divergence of these genes occurred long after the split of 
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N. munroi and the lineage leading to P. muelleri (Fig. 2), exclud- 
ing paralogy as an explanation for this grouping. Therefore, the 
transmission of these C 4 -optimized genes did not follow the 
species genealogy. A lateral transfer might have occurred via 
allopolyploidization, but this is unlikely given that none of the 
other nuclear markers isolated from P. muelleri showed affinities 
with those of N. munroi. It is more likely that the transmission 
occurred via hybridization, followed by backcrossing and intro- 
gression into the population, or alternatively through a lateral 
gene transfer event, as shown for the grass Alloteropsis (Christin 
et al., 20126), and suggested for the sedge Eleocharis (Besnard 
et al., 2009). Although the direction of the ppc-B2 gene transfer 
is not known yet, a transfer from N. munroi to P. muelleri is more 
consistent with the present data. In some regions of Australia, 
N. munroi and P. muelleri co-occur (Prendergast and Hattersley, 
1985; Supplementary Fig. S7), providing opportunities for a 
gene transfer that might have been favoured by natural selection. 
Some of the native P. muelleri ppc-B2 copies encode residues 
at positions 466 and 53 1 that generally characterize C 4 -specific 
PEPC (Besnard et al, 2009), which suggests that these genes 
previously fulfilled the C 4 function in P. muelleri and underwent 
a few changes toward C 4 optimization of the encoded enzyme. 
The acquisition of ppc-B2 genes from N. munroi that were 
more C 4 optimized would have been advantageous, leading to 
a replacement of the vertically acquired P. muelleri genes, and 
their subsequent pseudogenization (Fig. 3). 

Anatomical enablers of photosynthetic transitions in 
the Neurachninae 

The optimization of C 4 biochemistry is hypothesized to be one 
of the final steps in the evolution of C 4 plants (Sage et al., 2012), 
with the initial steps involving the reduction of the mesophyll 
cell to bundle sheath cell ratio, the enhancement of bundle 
sheath organelles, and the appearance of a photorespiratory C0 2 
pump, which characterizes nearly all known C 3 -C 4 intermedi- 
ates (Hylton et al, 1988; Griffiths, 1989; Hattersley and Watson, 
1992; Sage et al., 2012). The C 4 P. muelleri and N. munroi and 
the C 3 -C 4 N. minor have very similar foliar anatomies, with short 
interveinal distances, low mesophyll to bundle sheath tissue vol- 
ume, and high numbers of bundle sheath organelles, typical of C 4 
taxa (Hattersley et al, 1982, 1986; Brown and Hattersley, 1989). 
The C 3 taxa that separate them in the phylogeny have a compara- 
tively larger interveinal distance and smaller bundle sheath cells 
(Hattersley et al, 1982, 1986; Macfarlane, 2007), although the 
interveinal distance tends to be shorter than that of typical C 3 
taxa in N. queenslandica and N. annularis Macfarlane, and espe- 
cially inN. tenuifolia S.T.Blake, where additional minor bundles 
reduce the average interveinal distance (Hattersley et al, 1982; 
Macfarlane, 2007). C 3 species in the NeurachnelThyridolepis 
clade have conspicuous mestome sheaths, the tissue that has been 
recruited for carbon reduction in C 3 -C 4 and C 4 Neurachninae 
species (Dengler et al, 1985). Indeed, while this tissue is usually 
devoid of organelles in C 3 taxa, it contains some chloroplasts 
and a high number of mitochondria in these C 3 Neurachninae 
(Hattersley et al, 1982, 1986; Hattersley and Roksandic, 1983). 
The presence of these organelles in the mestome sheath may 
have facilitated the evolution of a photorespiratory C0 2 pump 



in the Neurachninae, enabling the subsequent transition to a 
complete C 4 pathway (Sage et al, 2012). As has been hypoth- 
esized for other anatomical characters in some eudicot lineages 
with multiple C 4 origins (Marshall et al, 2007; McKown and 
Dengler, 2007; Christin et al, 2011), the high organelle content 
of C 3 Neurachninae mestome sheaths may have contributed to 
the multiple origins of biochemical carbon-concentrating mecha- 
nisms in this grass group. 

Recurrent polyploidizations in the Neurachninae as 
diversifica tion facilita tors ? 

The distribution of ploidy levels in the species phylogeny (Fig. 1) 
suggests either multiple polyploidy events or a few followed by a 
return to a diploid state (diploidization). The latter suggestion can 
be disregarded because it would have led to a high chromosome 
number in the newly formed diploids, which is not observed in 
the Neurachninae (Table 1). In addition, multiple polyploidiza- 
tions are supported by the existence of intraspecific ploidy vari- 
ation in several taxa (Table 1; Prendergast and Hattersley, 1985). 
The data therefore suggest that seven independent polyploidiza- 
tion events have occurred in the Neurachne/Thyridolepis clade, 
one per polyploid species (Fig. 1). 

The analyses of markers from multiple genomes identified 
T. multiculmis and N. queenslandica as probable allotetraploids, 
the latter possibly representing a hybrid from a C 3 and a C 4 par- 
ent. In contrast, the tetraploid C 3 -C 4 intermediate N. minor does 
not appear to be the result of an allopolyploidization between 
C 3 and C 4 Neurachninae species, or at least not between known 
members, as the sequences from N. minor clustered together in 
the phylogenies at the same position for all plastid and nuclear 
markers (Fig. 1; Supplementary Figs S1-S4 at JXB online). With 
leaf anatomy similar to that of the C 4 N. munroi (Hattersley et al, 
1982), and a typically C 3 isotopic value (Table 1; Hattersley and 
Roksandic, 1983; Hattersley et al, 1986), consistent with min- 
imal or limited C 4 cycle activity (Hattersley and Stone, 1986; 
Hattersley et al, 1986; von Caemmerer, 1992), N. minor may 
represent a transitional stage toward an additional origin of full 
C 4 photosynthesis in the Neurachninae. As for N. minor, all the 
plastid and nuclear marker sequences belonging to the polyploid 
N. alopecuroidea, N. lanigera, and N. munroi clustered together 
in the phylogenies (Fig. 1 ; Supplementary Figs S 1-S6), suggest- 
ing that the dominant mechanism for increased chromosome 
numbers in this grass group was autopolyploidization. 

The two C 4 taxa and the C 3 -C 4 intermediate species of the 
Neurachninae are all polyploids of putatively recent origin 
(Fig. 1), and autopolyploidization might have allowed photo- 
synthetic diversification by providing additional copies of genes 
encoding enzymes that fulfilled housekeeping functions in the 
C 3 ancestors. A relaxation of selection due to genetic redun- 
dancy has indeed been suggested to facilitate the recruitment of 
enzymes into novel photosynthetic functions (Monson, 2003; 
Wang et al, 2009); however, as yet no empirical evidence exists 
to support this hypothesis, and comparison of distantly related 
grass genomes questioned the importance of gene duplication in 
the evolution of the C 4 syndrome (Williams et al, 2012). The 
apparent correlation of autopolyploidy and derived photosyn- 
thetic types in the Neurachninae offers an opportunity to test the 
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link between genetic redundancy and photosynthetic diversifica- 
tion within a sound phylogenetic framework. 

Conclusion 

The Neurachninae have long been viewed as a potential model for 
understanding the evolution of C 4 photosynthesis in the grasses 
due to natural photosynthetic variation in the group; however, 
evolutionary inferences were hampered by the lack of a phylo- 
genetic framework. Using multiple markers from the chloroplast 
and nuclear genomes, it has been shown that C 4 biochemistry has 
evolved more than once in this grass group, and may have been 
facilitated by anatomical enablers, a high frequency of autopoly- 
ploidization, and transfer of C 4 -optimized genes between related 
taxa. Molecular dating indicates that these events most probably 
occurred during the late Pliocene, an epoch of relatively low 
atmospheric C0 2 (Gerhart and Ward, 2010), when the Australian 
continent became drier (Bowler, 1976). The present-day distribu- 
tion of extant Neurachninae (Prendergast and Hattersley, 1985) 
indicates that the group appeared in the red sand soils and stony 
rises of semi-arid west-central Australia. These observations are 
consistent with recent models of C 4 evolution that hypothesize 
that climate deterioration in the past 5 million years selected for 
a series of traits that predisposed taxa to evolve novel photosyn- 
thetic mechanisms to deal with the combined effects of severe 
stress and low atmospheric C0 2 (Sage et al, 2012). 
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